NASA Contractor Report 191185 


3 o 7 3 £ 7 

20 f 


Analytical Skin Friction and Heat Transfer 
Formula for Compressible Internal Flows 


Lawrence J. De Chant 
Sverdrup Technology, Inc. 
Lewis Research Center Group 
Brook Park, Ohio 

and 

Marc J. Tattar 

Ohio Aerospace Institute 

Brook Park, Ohio 


January 1994 


Prepared for 

Lewis Research Center 

Under Contract NAS3-25266 


MSA 

National Aeronautics and 
Space Administration 


(NASA-CR-191185) ANALYTICAL SKIN 
FRICTION AND HEAT TRANSFER FORMULA 
FOR COMPRESSIDLE INTERNAL FLOWS 
Final Report (Sverdrup Technology) 

20 p 


03/02 


N9A-25I73 
Unc 1 as 
0207557 



ANALYTICAL SKIN FRICTION AND HEAT TRANSFER FORMULA 

FOR COMPRESSIBLE INTERNAL FLOWS 

Lawrence J. De Chant 
Sverdrup Technology, Inc. 

Lewis Research Center Group 
Brook Park, Ohio 44142 

and 

Marc J. Tartar* 

Ohio Aerospace Institute 
Brook Park, Ohio 44142 


SUMMARY 

An analytic, closed-form friction formula for turbulent, internal, compressible, fully developed flow has 
been derived by extending the incompressible law-of-the-wall relation to compressible cases. The model is 
capable of analyzing heat transfer as a function of constant surface temperatures and surface roughness as well 
as analyzing adiabatic conditions. The formula reduces to Prandtl’s law of friction for adiabatic, smooth, axi- 
symmetric flow. In addition, the formula reduces to the Colebrook equation for incompressible, adiabatic, 
axisymmetric flow with various roughnesses. Comparisons with available experiments show that the model 
averages roughly 12.5-percent error for adiabatic flow and 18.5-percent error for flow involving heat transfer. 


SYMBOLS 

A',B' "locally" defined variables (no formal association, used for convenience) 

B law-of-wall variable, 5.5 

C locally defined constant 

Cy coefficient of friction, 2r M /p u 2 

c p specific heat (constant pressure) 

c v specific heat (constant velocity) 

Cj,c 2 integration constants 

D pipe diameter 

E locally defined constant 

H total enthalpy, h + (1/2 )u 2 

h enthalpy 

k + roughness height; inner law variable 

M Mach number 
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Pr Prandtl number 

p pressure 

R universal gas constant; radius 

Re Reynolds number 

r recovery factor (Crocco-Busemann relation law), Pr 1/3 

T temperature 

u averaged velocity 

u flow velocity 

m + velocity; inner law variable 

1 D 

v* friction velocity, (x^/p^,) 

y physical wall variable 

y + wall variable; inner law variable 

Y ratio of specific heats, c f /c v 

K Von Karman constant, 0.4 

X Darcy friction factor 

K ID relative roughness 

p absolute or dynamic viscosity 

v kinematic viscosity 

p density 

x shear stress 

Subscripts: 

av average value 

w wall value 


INTRODUCTION 

A general, closed-form analytical solution describing a particular flow might be considered the "ultimate" 
mathematical fluid dynamics model because of its inherent flexibility and simplicity of use. Unfortunately, 
retaining adequate physical information in such a model is usually in conflict with our ability to reduce the 
solution to closed form. Typically, we are forced to discretize the problem and transform it to a large system of 
equations that are amenable to numerical solution techniques. The rapid growth in the availability of high- 
performance computing platforms has permitted the numerical solution of these large systems of governing 
differential equations by computational fluid dynamics methodologies. However, semiempirical or integral 
models are still capable of providing closed-form engineering estimates of flow physics in a highly cost- 
effective manner. This efficiency is particularly desirable for preliminary design analyses, where a large design 
space may need to be characterized. The goal of the present study was to develop a closed-form friction formula 
for turbulent, internal, compressible axisymmetric or two-dimensional duct flow. These flows are of interest both 
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in their own right and as a simple limiting approximation for more complex internal inlet and nozzle flows 
(ref. 1). 

Applications include modeling of acoustic liners for high-speed civil transport (HSCT) nozzles and national 
aerospace plane (NASP) propulsion system flow-path modeling. Further, this type of analysis may be used to 
provide first-principle estimates of engine component efficiencies (losses) for cycle analysis tools, such as the 
NNEP (the NASA/NAVY Engine Program, ref. 2). In addition, because heat transfer analysis is available, this 
model may provide significant information for use in structural analysis problems. 

As indicated, the friction model assumes fully developed flow and can model surface roughness as well as 
constant-wall-temperature heat transfer. The model is relatively simple because it essentially extends the incom- 
pressible law-of-the-wall relation to include compressibility effects. The Crocco-Busemann relation, an approxi- 
mate solution to the energy equation, is used to model the heat transfer. Consequently, the Prandtl number is 
assumed to be approximately 1 for this analysis. By requiring the governing equations to reduce to the incom- 
pressible limit, necessary constants of integration can be evaluated, thus permitting a closed-form friction 
coefficient relationship to be determined. In order to gain insight into the model’s reliability, the results of the 
analysis are compared with available experimental data. 


ANALYSIS 

The basic methodology followed in our development involves two steps. First, a compressible extension to 
the Prandtl mixing-length hypothesis is integrated to yield an effective velocity relation. Temperature/density 
dependence is modeled by an approximate energy integral, the Crocco-Busemann relation, and the ideal-gas 
relation. Second, this effective velocity is averaged and applied to the fully developed flow relation, yielding an 
implicit equation for the fully developed skin friction coefficient. 

To begin our analysis, we consider the shear stress closure by using Prandtl’ s mixing-length hypothesis 
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The energy equation is approximately integrated by using the Crocco-Busemann relation (ref. 3) 


H = CjM + c 2 


( 1 ) 


( 2 ) 


assuming that Pr ~ 1 and the pressure gradient is small. Therefore, with the definition of total enthalpy 


H = h + 



h + 


u 

= c , u + c 0 

2 1 2 


(3) 


and assuming that 


h = 




(4) 
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equation (3) then becomes 


c P r + — = c l» + c 2 


(5) 


The "no slip" boundary condition dictates that n(0) = 0. Evaluating equation (5) at the wall and assuming 
that the velocity u approaches an average as the temperature T goes to its average, the constants are computed as 


c, = 


c 2 = c p T w 


c p (T , ay - TJ + uj 2 


u. 


av 


( 6 ) 


Substituting into equation (5) gives 


~ . u 
c n T + — 
p 2 




v + -£h=- 


+ c p T w 


(7) 


And solving for T yields 


T = 


u 


av 


T — T -f 
av w 2c 


+ T -*L±- 
w 2 c„ 


( 8 ) 


Note that to partially ameliorate the Prandtl number (Pr = 1) assumption, a recovery factor r was added to give 
somewhat better accuracy. Its value depends on the local Prandtl number and for air in turbulent flow is about 
r = 0.88 (ref. 2). 


With the calorically perfect gas assumption 


and 





RTy 


(9) 

( 10 ) 


it can be shown that 
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Multiplying equation (11) by T iV /T w yields 


- 1 


LV 


+ - v-f-K 

L 1 W 


u 2 

av 


1 nw 2 ^av 

^r(Y - — 


( V 

u 


r + 1 




( 12 ) 


Applying the equation of state (ideal gas) and relevant boundary layer assumptions (dp/dy small) 


dp „ nPT Pw_ T 

-X- " °> Pav = /V p = P^ r ’ ~7T ~ ~T~ 

<ty 9 l w 


(13) 


combined with 
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yields 
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This expression is substituted into the original mixing-length formula 
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which yields 
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Separating and integrating 
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yields 


v* , M av 

In y + const = — _arcsin 

k # 


-2 B' 2 — + A' 

u av 


(45^ + A) 


/ 2 x l/2 


(19) 


where friction velocity is defined as 


v* 


f „ M/2 

T 


P W 


( 20 ) 


To evaluate the integration constant, the constraint applied is that the preceding equation must reduce to the 
incompressible case 


« + = J. In y + + B 

K 


( 21 ) 


for low Mach number flow (M av = 0). The integration constant may be split into two pieces for clarity: 


const = const j + const 2 


( 22 ) 


Rewriting the left-hand-side length scale (which is justifiable because it involves merely adding a constant) 
gives 


+ yv* 

>' y = it 

w 


(23) 


and from the incompressible law of the wall, const j is defined as 


const j = 


IlnIL 
K V„ 


= B 


(24) 
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The second constant const 2 is defined by considering the limit 5' defined previously as approaching zero. The 
term 


-2B' 


2 « 


lim^^ arcsin 

B' (4 5' 2 + A' 2 ) 1/2 


is bounded, whereas 


li m g'_ ) 0 


arcsin 
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(25) 


(26) 


is not. Since physically the limit should be finite, it is justifiable (because, again, it is merely a constant) to 
define the second integration constant as 


constj 


1 

— arcsin 
B' 


A' 

(45' 2 + A' 2 ) m 


(27) 


This rather subtle but physically reasonable line of argument defines the integration constants. 

Further, a classical modification of the incompressible law of the wall may be used to develop an expres- 
sion that takes the duct roughness into account. To do so. White (ref. 3) suggests the following modification to 
the incompressible law of the wall: 

( 

a* = I b. ZH + B 

K V„, 


where 


5 = 5.5 - — In (1 + 0.3 fc + ) 

K 



(28) 


Combining the previous relations yields the Van Driest effective velocity (ref. 2) 


v* 1 
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Because all flow properties must be related to either wall values or flow average quantities, the average of 
the preceding relation must be obtained. By averaging equation (29) over the area, it can be demonstrated that 
for a pipe the left-hand side becomes 

T = 1 In — + 1.75 - 2.5 In (1 + 0.3*+) (30) 

K v w 


The last term in this equation is the roughness term. Note that for smooth walls (k + = 0) this roughness term 
vanishes. 

Unfortunately, no such "elementary" integration will be available for the right-hand side (RHS). However, a 
simple approximation is available that calls for defining this average by replacing u by « av , yielding 


RHS 
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Separately, skin friction is defined as 


(31) 


V - T c /P 
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av av 
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where Cj is the desired 
and (19) as 


coefficient of friction. Thus, the friction velocity can be expressed from equations (32) 


( \l/2 


= V* = 


(2) 1/2 c 1/2 
5 f 


f M/2 
Pav 


\ Pw J 


u 


(33) 


and from equation (13) 


Then we obtain 
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permitting us to write the Reynolds number relation 
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where R here is the characteristic duct radius or hydraulic radius. Evaluating viscosity by using Sutherland s 
relation 


Pav 

Tv 


1.505 


1 + 0.505. 


where 
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and combining with equation (36) give 
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Recalling that 
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then it follows that 


k + = Re 


k (2) 1/2 
~D~ 2~ 



f ^ 


T 

x av 

1.505 

T 

T 

W 

1 + 0.505 _H 


T 


*av 


(40) 


9 


Combining the previous relations yields 
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By substituting equations (38) and (40) into equation (41) and defining the following constants (for both 
axisymmetric and two-dimensional flows) 
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the final form of the equation is obtained: 


, 1/2 


i(T - 1) ulr 


U/2 


= const + 1.77 In {cj 12 Re) + 1.77 In £ - 1.77 In |l + 0.2121 R 


(43) 


where const = -0.6005 for axisymmetric flow or 1.5086 for two-dimensional duct flow. 


RESULTS AND DISCUSSION 

Equation (38) represents the final closed-form solution for the skin friction coefficient. It is immediately 
obvious that the equation simplifies the solution for isothermal systems where T w = f av . For smooth pipes or 
ducts the roughness term disappears. An initial comparison with Prandtl’s law of friction for smooth pipes 

1 = 2.0 log [Re(X) 1/2 ] - 0.8 

a) m 

shows that the model is consistent to within 2 percent over a range from Re = 10 3 to Re = 10 7 (see table I, from 
ref. 41). This consistency is really no surprise in that the model was developed to recover the incompressible 
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limit as described by Prandtl’s law. However, Prandtl’s law does not consider roughness, heat transfer, or com- 
pressibility effects. 

In adiabatic systems where 


Ll = 1 + I(Y - 1)M 2 

r„ 2 

it is more difficult to see how the formula reduces. A much simpler approach would be to numerically compare 
the closed solution with experimental data points as a function of TJT^. 

Experimental data points have thus also been compared directly with the model. The more diverse these 
experimental flow conditions are, the more fully the model can be tested. Therefore, the experiments of interest 
range from two-dimensional duct flow to axisymmetric flow, both subsonic and supersonic flows as well as 
adiabatic and nonadiabatic flows. Discretion must be exercised when comparing actual data points with this 
model. First, extracting data indirectly off graphs and charts may often lead to inaccurate data readings. Second, 
experimental data points obtained with the use of Preston tubes are generally regarded as having an uncertainty 
of 10 to 20 percent (ref. 5). Laser interferometer skin friction (LISF) techniques are generally regarded as 
having better accuracy (see ref. 5). 

Table II shows skin friction estimates predicted by our model and measurements made in a square duct 
(ref. 6). This comparison is made for adiabatic, smooth, two-dimensional duct flow. In addition, the data point 
represents an average experimental friction coefficient across the width of the duct. The comparison is encourag- 
ing, although if roughness were included, the results would probably be somewhat better. 

To show this sensitivity to roughness, table III tabulates skin friction coefficients as a function of roughness 
for an adiabatic flow. Nevertheless, with roughness effects neglected, the model generally predicts skin friction 
within errors of less than 20 percent. 

In figure 1 the data points are much more scattered. This scattering may be related to the facts that only the 
inlet conditions were given for this experiment. Local conditions were estimated by using an adiabatic, viscous 
loss analysis (ref. 1). The results were then fed into the friction model. Therefore, the data points are only esti- 
mates of the relevant properties. 

Figure 2 compares the model with friction measurements made in a U-tube (ref. 7) in a fully developed 
region. Less error is associated with this experiment, possibly because data were obtained by using laser inter- 
ferometry rather than traditional Preston tubes. Because of this, the model agrees within an assumed 10-percent 
margin of uncertainty. This comparison shows that the model reduces accordingly for low Mach numbers to 
predict the skin friction for incompressible flows (as expected). Further experiments run on a square duct reveal 
similar results (ref. 8). 

Supersonic, adiabatic pipe flow data (ref. 9) were also compared with the present model. Table IV shows 
data for four Mach numbers and the corresponding predicted values. Agreement is rather good for all cases. 

To fully test the model, not only friction data are needed. Wall and flow temperatures must be known to 
compare the model’s heat transfer capability. Fortunately, several experiments have been done in this area. The 
first is shown in figure 3 (from ref. 10). Note that the 7' av /7' H , ratio of 1 indicates not an adiabatic system but an 
isothermal system for that particular location. Table V illustrates the model’s capability of predicting flows with 
roughness. 
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The recovery factor (r = 0.88) has also been used in the Crocco-Busemann relation. The value of the 
roughness term was chosen rather arbitrarily because no roughness figures were explicitly given. It was 
estimated to be k/D = 0.000095. However, White (ref. 3) suggests that the value chosen represents the point 
where roughness begins to play a significant role in heat transfer and may be considered a minimum value. 
However, other values of roughness may be more appropriate. Finally, the average Mach number was sonic 
(M = 1). 

With the ability to model surface roughness, the model may be further compared. By solving the friction 
formula for adiabatic, incompressible flows for various roughnesses, the model should reduce to the Colebrook 
equation (Moody’s chart). Figure 4 compares four roughness plots of the model with the Colebrook equation. 

As expected, the model recovered the Colebrook equation for small Mach numbers. In spite of this relatively 
good comparison, it is worth noting that the Moody chart is accurate to ±15 percent for design calculations 
(ref. 11). Figure 5 is included to show how the model extends the Moody chart to compressible flow for various 
Mach numbers. 

It is worth noting that figure 5 clearly presents the general trend that for fixed Reynolds numbers increas- 
ing Mach number tends to decrease skin friction coefficient values. Compressibility effects on the internal skin 
friction coefficient are directly presented in figure 6. Although it is not the intent of this report to discuss the 
detailed physics of this trend in any great detail, it is desirable to summarize the basic phenomenon. Schlichting 
(ref. 4) notes that for adiabatic, compressible flow the density must decrease strongly near the wall (by constant 
pressure and state), thus causing considerable boundary layer thickening. This reduction in the local gradient 
apparently reduces the skin friction coefficient. Unfortunately, our understanding of turbulent boundary layers is 
still rather phenomenalistic and can offer only limited insight into these problems. 

Solving for Cj in equation (39) may easily be performed by using various iteration (fixed point) techniques. 
A short routine using Newton’s method was used in this study. 

Overall comparison with experimental data yields no observable general trends of the model that may be 
identified as sources of error. It is worth noting though that because the model is based on Prandtl’s turbulence 
model, the present analysis cannot model the thin laminar sublayer (in fact, it is completely neglected). There is 
evidence that at low Reynolds numbers this effect becomes significant (ref. 3). White has demonstrated that for 
incompressible flow at Reynolds numbers of approximately 4000 the effect of the laminar sublayer becomes 
significant. 


SUMMARY OF RESULTS 

An analytic skin friction model for compressible, turbulent, internal, fully developed flow involving heat 
transfer has been developed by extending the incompressible law-of-the-wall relation to compressible cases. The 
formula recovers Prandtl’s incompressible law of friction for pipes (within 2 percent) for incompressible flow. 

In addition, the model shows good correlation with the Moody chart for similarly low Mach numbers. The model 
also shows good correlation with available data relating heat transfer and skin friction. Assuming the conditions 
under which the Reynolds analogy is valid, the skin friction can be directly related to the heat transfer. 


RECOMMENDATIONS 

In closing, it is worth noting that for these high-speed experiments the local properties (fluid density, 
Reynolds number, air speed) are difficult to measure. Often, they can only be estimated. As a result, agreement 
with available experiments is encouraging, but the experiments themselves may introduce some uncertainty. 
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Further, in contrast to the experimental work done in the area of flat plates, far fewer internal flow experiments 
have been performed to fully test this model. It is hoped that additional experimentation as well as agreement 
with numerically implemented turbulence models will prove this model to be accurate for rapid engineering 
estimation purposes despite its simplicity. 
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TABLE L— MODEL PREDICTION VERSUS PRANDTL’S LAW 


(INCOMPRESSIBLE) 


Reynolds 

number. 

Re 

■ 

Relative error. 
Mode! - Prandtl’s law 
Prandtl’s law, 
percent 

Model, 
X = 4C } 

Prandtl’s 

law, 

X = 4Cjr 

5.0x1 0 4 

0.0208 

0.0209 

0.48 

l.OxlO 6 

.0115 

.0116 

.86 

2.0x1 0 6 

.0102 

.0104 

1.92 

l.OxlO 7 

.0080 

.0081 

1.23 


TABLE H — MODEL PREDICTION VERSUS EXPERIMENTALLY 
MEASURED (ref. 6) SKIN FRICTION COEFFICIENT FOR 
SMOOTH, ADIABATIC, TWO-DIMENSIONAL 


DUCT FLOW 


Reynolds 

number, 

Re 

1 

Relative error, 
percent 

Mode! 

Experiment 

9.00x10 s 

l.SOxlO' 3 

1.80xl0 -3 

16 

l.OOxlO 6 

1.30xl0' 3 

1.30xl0" 3 

0 

2.00x1 0 7 

8-OOxlO -4 

l.OOxlO -3 

20 


TABLE m — MODEL PREDICTION 
FOR VARIOUS ROUGHNESSES 
IN TWO-DIMENSIONAL DUCT 


[Reynolds number, l.OxlO 6 ; 
Mach number, 3.O.] 


Skin friction 

Relative 

coefficient. 

roughness. 

c / 

K ID 

0.00190 

0.001 

.00160 

.0001 

.00158 

.00005 

.00156 

.00001 
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TABLE IV.— MODEL PREDICTION VERSUS EXPERIMENTALLY 
MEASURED (ref. 9) SKIN FRICTION COEFFICIENT 


FOR ADIABATIC PIPE FLOW 


Reynolds 

number, 

Re 

Mach 

number, 

M 

Skin friction coefficient, 

c / 

Relative 

error, 

percent 

Model 

Experiment 

8.00x10 s 

2.06 

2.75X10" 3 

2.80xl0 -3 

■9 

3.80X10 4 

2.84 

4.20x1 0 -3 

3.95xl0 -3 


4.50x10 s 

3.14 

2.00xl0 -3 

2.45xl0 -3 


2.25x10 s 

3.87 

2.30xl0 -3 

2.60x1 0 -3 

11.5 


TABLE V — MODEL PREDICTION VERSUS EXPER- 
IMENTALLY MEASURED (ref. 10) SKIN 
FRICTION FOR ISOTHERMAL PIPE 


FLOW WITH ROUGHNESS 


Temperature 

ratio, 

Reynolds 

number. 

Re 

Skin friction 
coefficient 
ratio, 

^/,prw/^fjneas 

0.909 


1.227 

1.00 


.960 

.476 

■ 

.875 

.400 


1.030 
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Skin friction coefficient, C, ^ Skin friotion coefficient, C, 


Mach 

number, 



106 10 7 108 


Reynolds number, Re 

gure 1 . — Skin friction coefficient versus Reynolds number for 
various Mach numbers. 


Temperature 

ratio, 

^av^vv 



Modified Reynolds number, Re 

Figure 3. — Skin friction coefficient versus Reynolds number for 
smooth pipe flow with heat transfer, Mach number, 0.6. (From 
reference 1 0.) 



Figure 2. — Skin friction coefficient versus Reynolds number for 
U-tube. Mach number, 0.1 . 


Relative 

roughness, 

kID 



0.001 


.0001 


.00001 


Smooth 


Modified Reynolds number, Re 


Figure 4. — Comparison with Moody chart. 
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iction coefficient, Cf 








Relative 

roughness, 

kID 


0.001 

.0001 
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